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1 Introduction 

A relationship between Painleve systems and infinite-dimensional integrable 
hierarchies has been studied. In a recent work jU [13] , we derive a coupled 
Painleve VI system, which admits a hypergeometric solution [bH [T6]. from 
the DS hierarchy of type A. In this article, we give a g-difference analogue 
of the above result. Namely, we propose a g-analogue of the DS hierarchy of 
type A and derive a g-Painleve system which admits a particular solution in 
terms of the g-hypergeometric function. 

The DS hierarchies are extensions of the KdV hierarchy for affine Lie 
algebras [2]. For type A^_ 1 among them, the hierarchies are characterized 
by partitions of a natural number N [5]. In this article, we propose a g-DS 
hierarchy corresponding to a partition (n m ) of mn. Note that its formulation 
is based on two preceding works [11], equivalent to the case m = 1, and [8], 
equivalent to the case n — 1. 

The coupled Painleve VI system given in [13] (or [15]) is derived from the 
DS hierarchy corresponding to a partition (n, n) by a similarity reduction. 
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Therefore, in this article, we investigate a similarity reduction of the g-DS 
hierarchy corresponding to a partition (n, n); we denote it by q-P( n ^ n y The 
system q-P( 2)2 ) coincides with the g-Painleve VI equation given by Jimbo and 
Sakai |7J . Thus we can regard q-P( n>n ) as a higher order generalization of the 
g-Painleve VI equation. The system g-P( n n ) also admits a particular solution 
in terms of the g-hypergeometric function n (f> n -i. 

The g-Painleve VI equation is a system of g-difference equations 



x(t)x(q-H) (y(q-H) - t^){y{q-H) - t0 2 

y{t)y{q- l t) 



0z0i 

where the parameters satisfy 



(x(t) — tc£i)(x(t) — ta 2 ) 
(x(t) — a 3 )(x(t) — a 4 ) ' 



PA 



q 



3 4 a 3 a 4 

It is given as the compatibility condition of a system of linear g-difference 
equations 

V(g _1 z, t) = A(z, t)Y(z, t), A{z, t) = Ao(t) + Ai(t)z + A 2 (t)z 2 , 
z(zI + B (t)) 



Y(z,q-H) 



-Y(z,t). 



;i-2) 



[z — q 1 ta\){z — q Ha 2 ) 
The coefficient matrices Ai(t) (i = 0, 1, 2) satisfy 



A 2 (t) 



Ki 
K 2 



Ao(t) has eigenvalues t9i,t9 2 



det A(z, t) = K\K 2 (z — tai)(z — ta 2 )(z — a 3 ){z — 0:4), 



where 



q 



«2 



a\a 2 



a.\a 2 



01 02 

The dependent variables x(t) and y(t) are given by 

(A(i))i2 g(x(t) - - to 2 ) 



x(t) = - 



(Ai(t)) 



12 



y(t) 



A(x(t),t) 



We derive the system ( II. ip from q-P( 2 ,2) with the aid of a g-Laplace transfor- 
mation for a system of linear g-difference equations. 

The g-hypergeometric function n </>n-i is defined by the power series 



nYn-1 



Oil, 

01,..., n -l 



q,t 



(oi; g)fc ■ ■ ■ («n-i; g)fc(«n; g)^ 



k=0 



{0i,q)k-. ■ (0n-i,q)k(q;q)k 
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where (a;q)k stands for the g-shifted factorial 

(a;g)o = l, (a; q) k = (1 - a)(l - qa) . . . (1 - g* -1 a) (k > 1). 



Note that it reduces to the generalized hypergeometric function n F n _i via 
a continuous limit q — > 1. We see that = „0„_i satisfies a n-th order 
linear g-difference equation 



where T q>t stands for a g-shift operator such that T qtt x(t) = x(qt). 

This article is organized as follows. In Section [2l we formulate a g-DS 
hierarchy and its similarity reduction corresponding to a partition (n m ) of 
mn. In Section [3J an explicit formula of g-P( n n ) is presented. We also 
discuss a group of symmetries and a Lax form for g-P( n n ). In Section HI the 
g-Painleve VI equation is derived from g-P(2,2)- in Section we show that 
q-P(n,n) admits a particular solution in terms of n (fin-i- 

Remark 1.1. The DS hierarchy corresponding to a partition (n m ) have been 
studied in [5] from a viewpoint of a Hamiltonian structure. It seems to be 
reduced from the q-DS hierarchy given in this article via a continuous limit. 

Remark 1.2 Q17J). The q-Painleve VI equation (11.11) is also derived from 
q-UC hierarchy. 

2 g-DS hierarchy 

In this section, we formulate a g-DS hierarchy and its similarity reduction 
corresponding to a partition (n m ) of mn. 

Let A and T\ be a grading parameter and a corresponding g-shift operator 
such that T\(\) = q\. We define mn x mn matrices ej, fj, hj (j — 1, . . . , mn) 



where E it j stands for a matrix with 1 on the (i,j)-th entry and zeros else- 
where. We also set 



((1 - q- 1 /3 1 T q , t ) ... (1 - g- 1 /3 n _ 1 T g ,i)(l - T g , t ) 

- t{\ - aiT q .t) ... (1 - a n _iT gjt )(l - a n T q Mx(t) = 0, 



by 




(j i mZ), 

(j G mZ; j mn), 



e j,k — ZjZj+l ■ ■ ■ Zj+k-l, fj,k — fj+k-i ■ ■ ■ fj+lf; 
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where ej +mn = ej and fj +mn = fj- Note that 

T\(e jtm ) = qe j>m , T A (/j, m ) = q~ X fj, m (.7 = 1,..., ran). 

Let us consider a set of grade one matrices which are mutually commutative. 
Explicitly, we take matrices 

n 

A; = 2J e i+C?-i)™>™ (i = 1, . . . , ra). 

5=1 

Then we have 

AjAj = AjAj = O, T x (Ai) = qAi (i,j = l,...,m;i^ j). 

Note that such choice of matrices is suggested by the previous work [TJ [T3J |5] ■ 
In the following, we formulate a g-DS hierarchy by using such matrices. 

Let ti and Tj (i — 1, . . . , m) be independent variables and corresponding 
g-shift operators such that 

Tifa) = qt h Ti(tj) = tj (j ^ i). 

Under a condition \q\ > 1, we consider a g-Sato equation 

T^Z- 1 = T l {W)(I-et l A l )W- x (i = l m), (2.1) 

with matrices of dependent variables 



van 00 van oo 



j=i i=i fc=i j=i fc=i 

where J stands for the identity matrix and e = 1 — q. We assume that 

^1,0^2,0 • • • z mn,0 = 1- 

(2.2) 

Under the system (I2.ip . we set 



* = W]\l[{I-q- k - 1 et i A i ). 



i=l fc=0 

Then we obtain a system of linear g-difference equations called a Lax form 

T i (*)=B i * (i = l,...,ra), (2.3) 

where 



The compatibility condition of the system (I2.3P is equivalent to 

Ti{Bj)Bi = Tj(Bi)Bj {i,j = l,...,m). (2.4) 

We call a system of g-difference equations (12. 4p a g-DS hierarchy. 

Next, we formulate a similarity reduction of the g-DS hierarchy. Under 
the system (12. 4p . we consider the following equation: 

m n 

T X (W) = q"T 1> ..., m (W)q- p , P = J2Y1 Mi+O'-Dm, (2.5) 

i=l 3=1 

where T^ m = T{T i+ i . . . T m and p\ + . . . + p m = 0. Note that 

pAi = Kip (i = 1, . . . ,m). 

We set 

tm oo \ 
i=l fc=0 / 

Then we obtain a Lax form 

T A (#) = M*, T0) = B^ (z = l,...,m), (2.6) 

where 

M = q p T2 j ... !m (Bi)T 3t ,,^ m (B 2 ) . . .T m {B m -i)B m . 
Note that the matrix M is also given by 

M = T x (W)q p (l - etiA x )(l - et 2 A 2 ) ... (1 - ^A^lf" 1 . 

The compatibility condition of the system (I2.6P is equivalent to 

Ti(M)Bi = T\(Bi)M, Ti(Bj)Bi = Tj(Bi)Bj = 1, . . . , m). (2.7) 

We call a system of g-difference equations (12.71) a similarity reduction of the 
g-DS hierarchy. 

Remark 2.1. The q-Sato equation (12.11) is derived as follows. Consider a 
matrix-valued function G defined by 

(m oo \ 
H H(I - q-^eUK) G(0, . . . , 0)[A, A" 1 ]. 
i=l fc=0 / 

Then we obtain a system of q- difference equations 

Ti(G) = (I - eUAi)G (z = l,...,m). 



It implies the system (12. ip via a triangle decomposition G = W~ 1 Z. In 
addition, a time evolution of the q-DS hierarchy is given in a framework of 
the function G by 

(m li oo \ 
i=l ji=l k=0 J 

Remark 2.2. Replacing M — >■ / — eM and Bi — > I — etiBi, we can rewrite 
the system (12.71) into 



[D x - M, A - Bi] = 0, [A - Bi, Dj - Bj] = (i,j = l,..., m), 
where 

A - : — , A - — t — • 

eX eti 

Via a continuous limit q — > 1, it reduces to the similarity reduction of the DS 
hierarchy of type A^ n _ Y given in 



3 Higher order g-Painleve system 

In this section, we present an explicit formula of the similarity reduction 
(12. 7p corresponding to a partition (n, n), where n > 2, in terms of dependent 
variables. We also discuss a group of symmetries and a Lax form for g-P( n n ). 
Consider a Lax form 

T A (^) = M^, T l (V) = B^ (i = l,2), (3.1) 

with 2n x 2n matrices 

M = T x (W)q p (l - e*iA0(l - et 2 A 2 )W-\ 
B i = T i (W){l-et i A i )W- 1 (2 = 1,2), 



where 



Ai = } j e 2j -i,2, A 2 = e 2 j,2, P = y^Pi(^2j-i - h 2j ), 
j=i j=i j=i 



and 

In oo 



j = l k=l 
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Then the compatibility condition of the system (13. ip is equivalent to the 
similarity reduction 

T 1 (B 2 )B 1 =T 2 (B 1 )B 2 , T i (M)B i = T x (B i )M (i = 1,2). (3.2) 

The matrices B±, B 2 and M can be expressed in terms of dependent variables 
Wj = Wj t \ and parameters Kj (j = 1, . . . , 2n). We assume that indices of Wj 
are congruent modulo 2n, namely Wj = Wj +2n . The equation (12.51) implies 

2i )2 (w 2 j-i) = q 2pl W2j~i, T 1>2 (w 2 j) = q~ 2px ~ X w 2 j (j = 1, . . . , n), 
from which we obtain 

2n n 

M = ^(1 - £Kj)hj + 2j(g pl etiw 2j - q~ px ~ x et 2 w 2 ,j_ 2 )t 2 j_x 

3=1 3=1 
n 

+ ^(q~ pl £t 2 w 2j+l - q Pl et x w 2 ^ x )e 2 j - q^et^ - q- pl et 2 A 2 , 

3=1 



and 



where 



2n 2n 

Bi = ^2 U i,3 h 3 + X] Vi >3 e 3 ~ £tiAi = X ' 2 )' 

3=1 3=1 



q ^(1 -enq-ij 



Ul,2j-1 = — ~ : 7^7 r, «l,2j = 1 + £tiT 1 (w 2 j-l)W2j, 

1 + g- 2 P J -- 1 et2W2j-2T 1 (w 2 j-i) 



ui,2j-i = etiw 2 j, v lt2j = -etiT 1 (w 2 j-i 

u 2 , 2 j-i = 1 + et 2 T 2 (w 2j - 2 )w 2j ^i, u 2 

V2,2j-i = -£t 2 T 2 {w 2j _ 2 ), V 2>2 j = £t 2 w 2j+1 
Note that the equation (12. 2p implies 



= g^(l -£K 2 j) 

2j 1 + g^et^-iT^; 



2n 2n 



Y[(1-6K J ) = 1, l[u hJ = l (2 = 1,2). 

We now assume that t 2 = 1. We also consider a change of variables 

p = q n , t = q 2np n\, z = q- n ^ pl - n+i y 2 (e\) n , 



a 



Xi (t) = g-^-^t^'+Vy-!, = g y - 1)(2 " 1+1 V l£ «; 2j , 

for j = 1, . . . , n. Replacing p — > q, we can describe the similarity reduction 
(13.21) as follows. 
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Theorem 3.1. The dependent variables Xj(t), Vj{t) (j — 1, . . . , n) satisfy a 
system of q- difference equations 



Xj(t) — Xj-i(t) - 



Vj{qt) - Vj-i{qt) 



1 + x^qtjyj^t) 1 + x i _ 1 (g*)y i _ 1 (i) ; 



(3.3) 



1 + Xjiqtfyjtt) 1 + Xjiqtfyj-tit) ' 
for j = 1, ... ,n, where 

bo = Q^bn, x (t) = tx n (t), y (t) = q~H~ l y n {t), 
with relations 

TT l + Xj(gt)y 3 -(t) = (n _i )/2 

We denote the g-Painleve system ( 13. 3p by g-P( n n ). As is seen in the next 
section, q-P( 2 , 2 ) coincides with the g-Painleve VI equation (II. ip . Note that 
q-P(n,n) reduces to the coupled Painleve VI system given in [131 EE] via a 
continuous limit q — y 1. 

The system g-P( n >n ) admits the affine Weyl group symmetry of type A 2 „_i- 
Let r,- (j = 0, . . . , 2n — 1) be birational transformations defined by 



bj-i ~ aj 



r 2 j-2(dj) = bj-i, r 2 j_ 2 (6 i _i) = a 5 , 

r 2j _ 2 (^_i(i)) = Zj-iW, r 2j _ 2 (yj-i(y)) = Vj-i(t) + ^ ^ _ ^ 
r 2 j- 2 (ai) = Oj, r 2i _ 2 (6j_!) = V^, 

r 2j - 2 (zi_i(t)) = r 2j -2{Vi-i{y)) = Vi-i{t) (i ^ j), 

and 

r 2i-i(«i) = bj r 2j -i(b j ) = aj, 

d ' — b m 

r 2j -i{x j (t)) = Xj(t) + J _ 3 r 2 j-i{yj{t)) = yj(t), 

yj\ l ) yj-i\ l ) 

r 2 j-x{ai) = a u r 2j - 1 {b i ) = b h 

r 2j -i(xi(t)) = Xi(t), r 2j -i(yi(y)) = yi(t) (i ^ j). 

Then q-P^ n , n ) is invariant under actions of them. Furthermore the group 
of symmetries (r , . . . ,r 2n -i) is isomorphic to the affine Weyl group of type 

/1 2n-l" 
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Consider a gauge transformation for the Lax form (13.11) 

n 

3=1 

Then we obtain a Lax form for q-Pt njV \ 

*(qz, t) = M(z, 0(z, t), *(z, qt) = B(z, 0(z, t), (3.4) 
with 2n x 2n matrices 



n 



M(z,t) = J2 a J E 2j-i,2j-i + J2 b i E W + (^(*) ~ ?~ lr V(*))#i,2 

3=1 3=1 

n—l n 

+ 2(^+1 (*) - ^(f^-E^y+i + ^2(Vj(t) - y j -i{t))E 2j - 1 ,2j 

3=1 3=2 

2n-2 

+ (Xi(t) - tX n (t))zE 2n ,l ~ 22 E 3-j+ 2 ~ tzE 1n-\.\ ~ zE 2n ,2, 

3=1 

and 



n n—l 

• £//.;(/ )/•-'•_>; . 2, - } j x j (qt)E 2 j t 2j+i - tx n (qt)zE 2 n,i 

3=1 3=1 
n-1 

— E 2 j-x j2 j + i — tzE 2n -i,i, 

3=1 

where i^j stands for the matrix unit. The group of symmetries defined above 
is derived from gauge transformations for the Lax form (J32 



rj (V(z, t)) = Ri(z, t)*(z, t) (j = 0, . . . , 2n - 1), 

where 

E> ( 4-\ T i — qCLi _j 

W*)- /+ Il(f )_ tlBW 2 ^ 

i^-- a (^ *) = / + 0- = 2, • • • , n), 

Note that a construction of those transformations is suggested by the previous 
works EDIIEI. 
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4 g-Painleve VI equation 

The system q-P( 2 ,2) is given as the compatibility condition of the Lax form 
$> 4 (qz, t) = M 4 (z, t)**(z, t), * 4 (z, qt) = B 4 (z, t)tf 4 (z, t), (4.1) 

where 



M 4 {z,t) 



V2(t) 
qt 



-1 



B 4 {z,t) 



at yi(t) - 

6i 
-fcz 

l + -xi(qt) 
-tz a — 



-tx2(qt)z 









x 2 (t) - xi(t) -1 

«2 y 2 (t)-y 1 (t) 
6 2 

-1 


l+X2{qt) V i(t) 2/2(0 

l + x 2 (g% 2 (t)J 



In this section, we derive the system (11.21) from the Lax form (14.11) with the 
aid of a g-Laplace transformation (cf. [SI IE])- 
First, we consider a gauge transformation 

%(z,t) = T 1 (z,t)S> 4 (z,t), 

where Ti(z,t) is a function such that 

n(qz,t) _ J_ Ti(z,qt) _ qt + x 1 (qt)y 2 (t) 
Ti(z,t) qa\ Tx(z,t) qta\ 

Then the Lax form (14. ip is transformed into 

*t(qz, t) = M* 4 (z, t)%(z, t), *l(z, qt) = B* 4 (z, t)%(z, t), (4.2) 

where 

M 4 (z,t) = M 4 (z,t), B 4 (z,t) = B 4 (z,t). 

qa± qta\ 

We set 

M^z, t) = M 4 * (t) + zM^t), B* 4 (z, t) = B* 4fl {t) + zB* 4>1 {t). 
Next, we consider a g-Laplace transformation 

$ 4 (C) - ^(rt) 



z%(z) -> 
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where e = 1 — q. Then the Lax form (14 .2p is transformed into 

$ 4 (g^C, t) = N A {(, t)$ 4 (C, *), $ 4 (C, 9*) = C 4 (C, *)$4& t). (4.3) 

where 

iV 4 (C,t) = (q-'eCI + M^t))- 1 (<M 4 * (t) + M^t)) , 
C A (C,t) = Bl fl (t) + e" 1 ^;^*) (/ - N 4 ((,t)) . 

Denoting C" 1 by z, we can restrict the Lax form (14. 3 p to the one with 3x3 
matrices 

* 3 (qz, t) = M 3 {z, t)V 3 (z, t), tt 3 (z, gf) = ,63(2;, t)tt 3 (z, £), (4.4) 
thanks to the following lemma. 

Lemma 4.1. For each of the matrices iV 4 (£,t) and C 4 (C,i) 7 the first column 
is equivalent to the fundamental vector *[1, 0, 0, 0] . 

Remark 4.2. TTie Lax form ( 14. 4 p coincides with the one for the q-Ql 3 hier- 
archy given in [8]. 

In a similar way, we can reduce the Lax form (I4.4p to the one with 2x2 
matrices. We consider a gauge transformation 

**(z,t)=T 2 (z,t)* 3 (z,t), 

where r 2 (z,t) is a function such that 

r 2 (qz,t) _ ai r 2 (z,qt) _ qta x 

r 2 (z,t) qbi r 2 (z,t) (qt + x 1 {qt)y 2 {t)){l + x 1 {qt)y 1 {t))' 

We also consider a g-Laplace transformation 

Then the Lax form (14. 4p is transformed into 

^(q-'C, t) = N 3 ((, t)$ 3 ((, t), $ 3 (C, qt) = C 3 ((, t)$ 3 (z, t). (4.5) 

Denoting e 2 a\biC, by z, we can restrict the Lax form (14. 5 p to the one with 
2x2 matrices 

^(q^z, t) = M 2 (z, t)* 2 (z, t), * 2 (z, qt) = B 2 (z, t)V 2 (z, t), (4.6) 
thanks to the following lemma. 

11 



Lemma 4.3. For each of the matrices N 3 ((,t) and C^Ct), the first column 
is equivalent to the fundamental vector *[1, 0, 0]. 

The matrices M 2 (z,t) and B 2 (z,t) are of the form 
», / + x _ M 2|0 (t) + zM 2il (t) + z 2 M 2 , 2 (t) B 2i0 (t) + *fl 2 ,i(*) 

M 2 {Z,t) — — — , tf 2 {Z,t) — — 



(z-t)(z-l) 



z-t 



where 



M 2)2 (t) = - 



a 2 y 2 {t)-yi{t) 
b 2 

1 



l+x 2 (qt) yi (t) 

l+x 2 (gt)i/ 2 (^ 



l + zi(gt)3/i(t) 

In order to derive the system (jl.2p . we consider a gauge transformation 



y/ ^ = (9 X tz u ,q ^ooiqziqfc \r 3 {t) 



o T 4 (t) 



V2(t)-yi(t) 
a2~b'2 

1 



* 2 (z,*)- 



where r 3 (z, i) and r 4 (z,t) are functions such that 



Tz{qz,t) r 3 (z,qt) (1 + xi(gt)j/i(t))(l + x 2 (qt)y 1 (t)) 



T 3 (z,t) 

T 4 (qz,t) 



r 3 (z,t) a 2 
T±(z,qt) 1 + x 1 (qt)y 1 (t) 



r 4 (z,t) T 4 (z,t) 1 + x 2 (qt)y 2 (t)' 

Then the Lax form (14. 6p is transformed into 

Y(q-\t) = M 2 (z,t)Y(z,t), Y(z,t) = B2{z > q ~ H) Y(z,q-H). (4.7) 
Here the coefficient matrices satisfy 



M 2 {z,t) = M 2fl {t) + zM 2>1 {t) + z 2 M- 



2,2, 



M- 



2,2 



a 2 
b 2 



M 2 (t) has eigenvalues ta\,tbi 



and 



detM 2 (z,t) = a 2 b 2 {z - t)(z - ai&ig~ 1/2 f)(z - l)(z - a^b^q 1 ' 2 ), 

B 2 (z,q- 1 t) = B 2fi (q-H) + zI 2} 

det B 2 (z, q~ l t) = (z - q~H)(z - q^aihq^H). 
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We can show them by direct computations; we do not state its detail here. 
By setting 

a x b x gV2 
«i = 1, a 2 = —rpr, a 3 = 1, a 4 = — — , 
q 1 '^ a 2 b 2 

Pi — TP?' ^ 2 = ~TF?' ^3 — — , At = -r, 
q L i z q L ' z a 2 b 2 

we arrive at 

Theorem 4.4. The matrix-valued function Y(z,t) solves the system (jl.2p . 

Corollary 4.5. Under the system q-P( 2 ,2), we set 

= t(x 2 (t) -Si(t))£i(t) , , = x 2 (gt)(qt + x 1 (qt)y 2 (t))i(}i(t) 
X[) Ut) ' V[) (l + x 2 (qt)y 2 (t))Mt) ' 

where 

£i(t) = qtx^y^t) - x 1 (t)y 2 {t) - QtB 2 (t)j/i(t) + x 2 {t)y 2 (t) - (6 a - ai)gt, 
6(t) = (*z 2 (*) - xi(t))(x 2 (t) - xi(t))G/ 2 (*) - g%(t)) 

+ (6i - a^qtxxit) + {(a 2 - - (a 2 - ai)}gix 2 (i), 
^(t) = (1 - a 1 6 1 g 1/2 t)x 2 (gt)y 2 (t) + gt - ai&ig 1/2 *, 
^ 2 (t) = o 2 (l - ai6ig 1/2 t)xi(g£):E 2 (gt)2/ 2 (t) 

+ ai(gt - a 2 biq 1/2 t)xi(qt) + (a 2 - a 1 )qtx 2 {qt). 

Then those variables satisfy the q-Painleve VI equation (II. ip . 

Note that g-P( 2j2 ) gives explicit formulas of x\(qt) and x 2 (qt) by 



xi(g£) 
x 2 (qt) 



+ (<Z* - axb iq VH) yi {t) - (1 - UihqVHfaity 

Z^W 

^4(t)2/i(t) + a 2 (gt - ax&xg 1 ^)^^) - a 2 (l - a^q^y^t) 

where 

6(0 = aig 1/2 t(x 2 (t) - xi(t))(y 2 (t) - yi(t)) - gt + aia 2 g 1/2 t, 
= (x 2 (t) - xi(t))(i/ 2 (t) - gfyi(t)) - 6i(gt - ai a 2 q l/2 t). 
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5 g-Hypergeometric function n (j) n -i 

In this section, we show that q-P( n .n) admits a particular solution in terms 
of the g-hypergeometric function n <f) n -x- 

Proposition 5.1. Under the system q-P(n, n ), we consider a specialization 

n 

Vj {t)=o (j=i,...,n), n%=? (n ~ i)/2 - 

3=1 

Then a vector of the variables x(t) = t [xi(t), . . . ,x n (t)} satisfies a system of 
linear q-difference equations 

x(g- 1 t)= (^o + TZ ^)x(t), (5.1) 

with n x n matrices 

n n n n n 

A o = Y h 3 E 3\3 Y Y [b J ~ "J )l ' : '-r A Y Y [a ' ~ h i) E hr 

j=l i=l j=i+l i=l j=l 

We always assume that 

a,j ^ Z, a,i — a,j Z (i, j — 1, . . . , n; i ^ j). 
Note that is an upper triangular matrix and 

n n i— 1 

In • .1: Y a > E ij + S Y { "' - h ^' : '-r 

3=1 i=l 3=1 

is a lower triangular matrix. 

We consider a formal power series of x(£) at t = 

oo 

x(t)=t 1 ^ a ^(g- 1 t) fc x fc . 

fc=0 

Recall that |g| > 1, namely < 1. Substituting it into the system (15. ip . 
we obtain 

(A + A x - a x I)xo = 0, (A + A x - a!g- fc /)x fc = A x fc _! (A; > 1). (5.2) 
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The matrices A + A\ — a\I and A + A\ — a\q~ k I are of rank n — 1 and 
n, respectively. Hence the recurrence formula (15.21) admits one parameter 
family of solutions. Its explicit formula is given by 



CL\ . . . Cb n 



(k>0), 



where 



3-1 



^0,1 — 1; ^0,, 



II (J = 2,...,»), 



1=1 



— a l 



and 



(gg; • • • , (gfc; g)*, (g; gk ■ ■ ■ , g)* 
?)*, • • • , <?)*, (^; g)*, ...,(£; g)*(g; g)* 

for > 1. Then we arrive at 

Theorem 5.2. The system ( 15. ip admits a solution 



-x 



j * 



x(t) = t 



log ai 



where 



3-1 

n 

i=l 



</?,(*) = n (f) 



n-1 



n.—L n 1 1 /y "'i 

y fo. ) • • • ) y a . , ) h - ) • • • 5 y fe n ? /> 

^ a 2 ' a j a j + l a n 



-i h 



ai . . . a n 



-g * 
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